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ABSTRACT 

Cosmic-ray interactions with interstellar gas produces both 6 Li, which accumulates in the interstellar 
medium (ISM), and ir° mesons, which decay to gamma-rays which propagate throughout the cosmos. 
Local 6 Li abundances and extragalactic gamma-rays thus have a common origin which tightly links 
them. We exploit this connection to use gamma-ray observations to infer the contribution to 6 Li 
nucleosynthesis by standard Galactic cosmic-ray (GCR) interactions with the ISM. Our calculation 
uses a carefully propagated cosmic-ray spectrum and accounts for 6 Li production from both fusion 
reactions (act — > 6 Li) as well as from spallation channels (p, a + CNO — > 6 Li). We find that although 
extreme assumptions yield a consistent picture, more realistic ones indicate that solar 6 Li cannot be 
produced by standard GCRs alone without overproducing the hadronic gamma rays. Implications for 
the primordial 6 Li production by decaying dark matter and cosmic rays from cosmological structure 
formation are discussed. Upcoming gamma-ray observations by GLAST will be crucial for determining 
the resolution of this problem. 

Subject headings: cosmic rays - gamma rays: theory - nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

Cosmic-ray nucleosynthes is is the only known Galac- 
tic source of the 6 Li l|Vangioni-Flam et alJ 119991 
IFields fc Qlivelll999h . Thus, it is a standard belief that 
the observed solar abundance of 6 Li was produced by 
Galactic cosmic-ray (GCR) interactions with the inter- 
stellar m edium (ISM), where aa — > 6 Li is the dominan t 
channel ijSteigman fc Walked 119921: iMontmerld Il977|) . 
However, hadronic CRs also produce gamma rays, and 
thus GCR inte ractions in normal galaxi es are guaranteed 
to contribute l|Pavlidou fe Fields! I2002|) to the observed 
extragalactic gamm a-ray background (hereafter EGRB; 
iStrong et al .1120 04). Moreover, since they both originate 
from CR interactions, 6 Li and hadronic gamma-rays are 
ti ghtly related. 

IFields fc Prodanovid 1)20051) established a simple and 
model-independent connection between lithium and "pi- 
onic" 7-ray production (pp — -> tt — > 77) by a given 
cosmic-ray population. Using this tool with simplifying 
assumptions gave the alarming result that the solar 6 Li 
abundance, if produced entirely by GCRs, demands a pi- 
onic 7-ray intensity exceedin g the entire observed EGRB 
ijFields fc Prodanovidl200 5r). Given the current interest 
in 6 Li, this result thus deserves a thorough investigation. 

In this paper we revisit the problem of Li-7-ray con- 
sistency with a more precise and realistic calculation. 
We now employ a carefully propagated CR spectrum, 
as opposed to the standard sing l e-pow er law spectrum 
adopted in IFields fc Prodanovid l)2005jh Moreover, in- 
stead of using a convenient fit fo r the pionic 7-ray spec- 
trum ( Pfrommer & Enfilin 2004) we now calculate it self- 
consistently from our CR spectrum. We also estimate the 
spallation p,a + CNO — ► 6 Li contribution to the solar 
6 Li abundance. These effects slightly reduce but do not 
eliminate the discrepancy. Moreover, the only remaining 
effect we expect to be important- 6 Li destruction as it 
is processed through stars-makes the problem more se- 
vere. The net effect is that in a realistic calculation, the 



observed EGRB allows for only w 60% 6 Li Q to be pro- 
duced by standard GCRs. Only a conspiracy of extreme 
assumptions gives GCR production of the solar 6 Li that 
does not at the same time saturate the observed EGRB. 

Our result represents a strong hint for the need of 
a new 6 Li source. Recent suggestions such as dark 
matter and low-energy cosmic rays are discussed in 
jpl Upcoming gamma-ra y observations by GLAST 
ijGehrels fc Michclson 1999) will better constrain (or de- 
termine!) the pionic 7-ray fraction of the EGRB and 
will thus be the key in determining the severity of this 
problem. 

2. LITHIUM-GAMMA-RAY CONNECTION 

In IFields fc Prodanovid l)2005l) we formally demon- 
strated and quantified the tight connection between CR 
lithium synthesis and hadronic 7-ray production. We 
showed that both observables are a measure of the time- 
integrated CR flux (fluence F). Specifically, the ratio 
of the "pionic" 7-ray intensity I lw (integrated over the 
entire energy spectrum) and 6 Li abundance (baryon or 
mole fraction 6 Li = Yq = ri6/"-baryon) produced in fusion 
reactions with the ISM can be expressed essentially as 
the ratio of their reaction rates 
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This factorizes into a product of nuclear, cosmological, 
and cosmic-ray parameters, and a "Copernican" factor 
-Favg/Fiviw- Cosmology enters via the comoving baryon 
number density rib = 2.52 x 10~ 7 cm~ 3 . The cosmic- 
ray and ISM helium abundances are taken to be y™ = 
y„ m = 0.1 (jji = ni/nu)- The flux- averaged pionic 7- 
ray production cross-section is a n = 2£ Q ,£ 7r cr 7r o where the 
factor of 2 counts the number of photons per pion decay, 
avo is the cross section, £„. is the pion multiplicity, and 
the factor = 1.45 accounts for pa and aa reactions 
(perme r 198G). For ct 6 we have used a recent result of 
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iMercer et al.l l)2001l) . which for 6 Li differs significantly 
from the old values. The use of the new cross sections 
results in a lower 6 Li production. 

The ratio -Favg/^MW is the ratio of the line-of-sight 
baryon-averaged fluence to the local fluence; this "Coper- 
nican ratio" compares the cumulative cosmic-ray activity 
of our Galaxy to that of an average star-forming galaxy. 
Following our previous work we will initially take this 
factor be ~ 1, i.e. that the Milky Way CR flux through 
out the history can be approximated with the cosmic 
mean. We will then examine the consequences that this 
ratio differs significantly from unity. 

3. COSMIC-RAY SPECTRUM 

In eq. Q), the Li-7-ray proportionality depends on the 
ratio of the mean cross sections u-JdQ. These must 
be properly averaged over th e GCR energy spectrum. 
In iFields fc Prodanovici 1)20051) we have adopted a stan- 
dard propagated cosmic-ray spectrum which is a single 
power-law in total energy with a spectral index a = 
2.75 over the entire relevant energy range. While this 
is a commonly-used rough approximation to the GCR 
spectrum, it becomes inaccurate at energies < 1 GeV, 
where ionization energy losses dominate over escape 
losses. Because aa — > 6 Li threshold energy is at ~ 10 
MeV/nucleon, while pp ->■ tt° threshold is at ~ 280 MeV, 
the Li-7 connection is particularly sensitive to GCR be- 
havior at very low energies. Thus in this paper we refine 
on the analysis presented in lFields fc Prodanovici ( 2005) 
by calculating and implementing a car efully propagated 
CR s pectrum for a leaky box model ijMeneeuzzi et al.l 
1971). 

The resulting CR spectrum calculated using the stan- 
dard leaky box model, assuming a standar d source spec- 
trum that is a power-law in momentum l)Gaisserlll990l 
e.g.,). This gives a CR flux ~ 4 times higher around 
aa — > 6 Li threshol d , com pared to the one used in 
IFields fc Prodanovici J2005) where a single power-law 
spectrum was assumed, while for energies > 100 MeV 
a single-power law spectrum is a good approximation. 

4. PIONIC GAMMA-RAY SPECTRUM 

iPfrommer fe Enfilml l)2004|) provide a useful 
pa rametrization of the pionic 7-ray spectrum used 
in IFields fc Prodanovici 1(20051) . However, here we 
numerically calculate the pionic 7-ray spectrum in full 
detail, by adop ting the isobar+scaling model as given in 
iDermerl (^986); the pionic spectrum we adopt uses the 
same cosmic-ray spectrum as the 6 Li production, and 
thus is self-consistent. 

In order to calculate 7 77r one needs to know the history 
of the CR sources and the targets. Both histories come 
from the cosmic star-formation rate . As described in 
detail in IFields fc Prodanovici l)2005l) we can obtain the 
GCR pionic 7-ray spectrum integrated over the history 
of the sources (equation 26 of lFields fc Prodanovic 2005, 
same parameter values used) . The cosmic star- formation 
rate alone fixes the shape of the pionic EGRB, but re- 
quires a normalization that physically connects the star 
formation rate to the cosmic-ray flux, and which normal- 
izes the present gas fraction in a typical galaxy. In order 
to place an upper limit to the pionic EGRB, we allow 
this normalization to vary freely to maximize the pionic 
7-ray flux consistent with present EGRB observations 
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Fig. 1. — In the upper panel of this figure, we plot the pionic spec- 
trum (dotted green line - maximized, dashed blue line - normalized 
to the Milky Way) , compared t o the observe d EGRB spectrum 
(solid line, fit to data); we use the Strong et al. (2004) data points, 
which are given in red crosses. The bottom panel represents the 
residual function, that is, log[(7_B 2 ) b s /(/-B 2 )7r] = iog(I \ 3s / 



({Fields fc Prodanovidl2005l iProdanovic fc Fields! I2004J) . 

This is presented in Fig. 1 as a do tted green lin e . The 
observed EGRB spectrum is that of lStrong et alJ 1(20041) 
and is plotted as red data poin ts, with a black solid line 
fit ((Fields fc Prodanovici I2005D . Finally, we find maxi- 
mal pionic 7-ray fraction to be 58% of the total observed 
EGRB. 

More realistically, we can use the Milky Way to deter- 
mine both the scaling between the star formation rate 
and the cosmic-ray flux, and the present-day gas frac- 
tion. We do this following iPavlidou fc Fields] l|2002|) . us- 
ing a present gas content of the Milky Way M gaSi MW ~ 
10 10 M Q , and a star formation rate 3.2 M Q /yr. The 
resulting 7-ray spectrum is presented in Fig. 1 as 
a blue dashed line. This corresponds to the pionic 
7-ray contribution expected from the normal galax- 
ies. In addition to this guaranteed component to 
the EGRB , unresolved blazars will a lso contribute sig- 
nifica ntly l|Stecker fc Salamonl Il996t IPavlidou fc Fields! 
2002), presumably comprising much or all of the remain- 
ing signal. 

Having determined a upper limit and a more realistic 
estimate to / 7jr one can find the corresponding Li abun- 
dance, via eq. This is our main goal, to which we 
now turn. 

5. UPPER LIMITS ON AND ESTIMATES OF 
GCR-PRODUCED 6 Li 

In this section we calculate limits to and estimates of 
the 6 Li produced by GCRs that are allowed by preset 
EGRB data. We present our results in steps of increasing 
realism. For now we retain the Copernican assumption 
that the Milky Way cosmic-ray fluence is typical of star- 
forming galaxies (-Fiiw/^mw = 1); we will revisit this 
assumption in the final section. 

1. By combining with the maximal pionic 7-ray 
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fractio n and procedure described in lFields fe Prodanovicl 
( 2005), we find the fraction of 6 Li abundance produced in 
aa -» 6 Li re action to be 6 Lw = 0.61 6 Li (( 6 Li/H) Q = 
1.53xlO~ 10 : IAnders fc Grevessell989ft . This corresponds 
to an extreme upper limit for all 6 Li produced by the 
GCR aa reaction. 

2. Though the aa reaction with the ISM is the dom- 
inant channel for 6 Li production, a non-negligible con- 
tribution, especially at higher metallicity, comes from 
the spallation reactions p, a + CNO — > 6 Li (both for- 
ward and inverse kinematics, that is fast heavy nuclei, 
are included) . If the fusion and CNO reaction rates were 
to be equal the required oxygen abundance should be 
(0/H) cq = 0.51 (0/H)q. This now sets the normaliza- 
tion and allows us to calculate the total 6 Li abundance 
produced from all channels under the extreme assump- 
tion that the ISM was at solar metallicity over the Galac- 
tic history. We find that 6 LiccR = 1-79 6 Li Q , which now 
represents the extreme upper limit for all 6 Li produced 
by GCRs. 

3. Because the cosmic-ray CNO abundance is a direct 
function of the Galactic supernova rate, a precise calcu- 
lation introduces a factor of 1/2, that is, instead of as- 
suming solar metallicity through out history one should 
use an average value of (0/H) oq = 0.5 (O/H) . This re- 
sults in the total allowed GCR-produced 6 Li abundance 
of 6 LiocR = 1-20 6 Li©, which is still consistent with the 
standard picture. 

4. So far we have been taking the maximal (Fig. 1, 
dotted green line) pionic 7-ray fraction as allowed by the 
present EGRB data 1 , where we have (without justifica- 
tion) ignored the normalization and just used the shape 
of our spectrum. However, it is unrealistic to assume 
that the entire emission is due to GCRs. Indeed, inde- 
pendent of the details of our galactic 7-ray estimate, it is 
clear that the EGRB must contain a large and perhaps 
dominant contribution from unresolved AGNs (blazars) 
and so the galactic signal must leave room for this and 
cannot saturate the observed level. An estimate of the 
normalized GCR pionic 7-ray component of the EGRB 
(Fig. 1, dashed blue line) yields a spectrum that is a 
factor of 2.1 lower than the maximized value. Thus, in 
this most honest case, we find 6 L1gcr = 0.57 6 Li which 
now falls short by about a factor of 2 from a standard 
picture of cosmic-ray 6 Li nucleosynthesis. 

5. For inverse CNO kinematics a non- negligible LiBeB 
production comes from t wo-step spallation r eactions, eg. 
O + H B + H -> 6 Li lIKneller et al.ll20f)l . For exam- 
ple, the production rate of 6 Li from two-step reactions 
of fast oxygen is ~ 40% of single-step fast oxygen spalla - 
tion reactions, for a fixed A = 10 g/cm 2 l|Knellerl [2006'). 
However, when two-step inverse CNO kinematics is taken 
into account, the overall increase is only slight and the 
result now becomes 

6 Li GCR = 0.59 6 Li (2) 

Even in the most extreme assumption that the two-step 
rates are equal to the single-step inverse CNO kinematic 
rates, the resulting 6 Li abundance would still be only 
63% of the solar. 

1 Determination of the EGRB relies on the subtraction of 
the Galactic Plane and is thus model-dependent. Although 
IKeshet et alj 120041) report only a limit to the EGRB this would 
only strengthen our result. 



6. Finally, one has to remember that the observed solar 
6 Li abundance is not the total lithium abundance pro- 
duced, due to astration, that is, the fact that some of the 
gas was already processed by stars. Due to very fragile 
nature of this isotope, 6 L1q is only the lower bound on 
the total 6 Li produced. For a rough estimate of the level 
of astration one can use deuterium. It is well established 
that Big Ba ng nucleosynthesis is the only important 
sourc e of D ijEpstein et alJ 119761: IProdanovic fc Fields! 
120031) and that D is easily destroyed in stars due to a sim- 
ilarly fragile nature. Thus by compa ring the solar neb- 
ula D abundance D prcso i = 2.1 x 10 -5 (|Oeiss fc Gloecklerl 
1998) with the abundance determined from 5 best quasar 
absor ption systems Dqso = 2.78 x 10~ 5 (|Cvburt et alJ 
we find that roughly ~ 25% of the gas has passed 
through stars. Thus 6 Li Q is about ~ 75% of 6 Li tot , 
and our calculated GCR 6 Li now becomes 6 LiocR ~ 
0.45 6 Li tot . 

6. DISCUSSION 

In this paper we have used the connection between 6 Li 
and pionic 7-rays produced in CR interactions, in order 
to calculate the total allowed 6 Li abundance that can be 
produced by GCRs. We have used a CR spectrum, care- 
fully propagated according to the leaky-box model, while 
the pionic 7-r ay spectrum was calculated based on the 
iDermerl (^986) model. A realistic, detailed calculation 
that includes 6 Li production from both fusion reaction 
with the ISM and spallation CNO channels (2-step in- 
verse kinematics also included), yields a 6 Li abundance 
that is only w 60% of the total 6 Li produced, if standard 
GCRs are the only relevant source. Correcting for astra- 
tion will result in even lower 6 LiccR abundance at the 
level of ~45% 6 Li tot - 

Our result either indicates the need for a new impor- 
tant source of 6 Li beyond standard GCR nucleosynthesis, 
or it points to a possible failure of the usual assumption 
that the average interstellar GCR flux tracks the instan- 
taneous star formation rate. We consider each possibility 
in turn. 

Additional sources of 6 Li are of considerable cur- 
rent interest, because of the re cent report of a 6 L i 
plateau in metal-poor halo stars (jAsplund et al. 2005). 
As with the familiar 7 Li Spite plateau, an analogous 
6 Li feature would suggest a pre-Galactic source of 6 Li. 
And indeed, recently two very different cosmological 
sources of 6 Li have been proposed: (1) production 
in the early universe, stimulated by supersymmetric 
dark m atter particle decays dur i ng big bang nucleosyn- 
thesis flDimonoulos et alT 119881: Kawasaki et all 120051: 
.Tedamzik et alJ 120051 lEllis et all 120051 iKusakabe et alJ 
2006); and (2) production during the virialization and 
baryonic accretion of large -scale structures, w hich gen- 
erates cosmological shocks ( Mini ati et al.ll200"fH) that can 
in tu rn accelerate a popul ation of cosmological cosmic 
ravs llSuzuki fc Inoudbool lBlasill2004l but see lPrantzosl 
(2006) for constraints). 

The 6 Li plateau is < 10% 6 Li Q , and thus whatever its 
source is, it will not be able to account for the factor > 2 
discrepancy between 6 Li Q and 6 LiocR we have found. 
However, the existence of the 6 Li plateau at the 10% level 
of the solar abundance for metallicities [Fe/H] < — 1, can 
be used as a constraint to any non-standard 6 Li source 
that is expected to account for the potentially missing 



4 



~ 40% of 6 Li Q . Moreover, 6 Li plateau would indicate 
that such a source would have to become important only 
at late times, and near-solar metallicities. 

We note that another additional source of 6 Li could 
come from a population of CRs having low energies 
(< 100 MeV). Such particles are excluded from the so- 
lar system and hence not directly constrained obser- 
vationally. A large flux of such particles, well above 
the extrapolated observed high-energy trends, could pro- 
duce large additional amounts of 6 Li but no pions and 
hence no pionic 7-rays. Indeed, recent ob servations of 
Hjj~ in molecular clouds l)McCall et al]l200ll) seem to de- 
mand a large low-energy CR flux in the neighborhood 
of these clouds. On the other hand, low-energy CRs 
widespread enough to participate significantly in LiBeB 
nucleosynthesis on Galacti c scales face strong constraints 
that come from energetics llRamatv fcjjing ejJeltejl999|l 
and f rom LiBeB abundance Tedios ~^Van^m^^\amet^l\ 
1998). These limits are evaded if solar 6 Li reflects a 
localized low-energy CR enhanceme nt, either due to a 
hypernova-like Type Ic supe rnova (|Fields et all 12001 
iNakamura fc Shigevamall2004|) . o r to solar CR produc - 
tion in the protosolar nebula fe.g.. lGounelle et aTl l2006'). 
In either case, the other 7 LiBeB isotopes will be produced 
and constrain the allowable 6 Li contribution. 

In this work we have assumed that the Milky Way 
CR fluence can be approximated with the cosmic mean. 
Therefore, our result might indicate that more 6 Li was 
produced than 7-rays would suggest, which would be the 
case if the Milky Way CR flux was at some time(s) a 
factor of ~ 2 (on average) higher than the typical CR 
flux in a normal galaxy. 

If indeed 6 Li points to enhanced CR activity, this in 
turn would point to anomalies in Milky Way star for- 
mation and/or CR properties. We have assumed that 
the CR fluence here is typical of the mean star-forming 
galaxy (i^Mw/^avg = 1)- If instead our Galaxy had a 
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more vigorous CR history, this could account for the 
difference. Also, we have in our most realistic assess- 
ment used the present local cosmic-ray/star- formation 
ratio $mw/^mw; if this departs from the cosmic mean, 
this too could account for the 6 Li discrepancy. Both of 
these solutions have implications for (and quantify) the 
utility of the Milky Way and Local Group as representa- 
tive cosmological samples of star forming galaxies. 

Upcoming gamma-ray experiments will go far to clarify 
the nature of Galactic and extragalactic pionic gamma- 
rays, and hence 6 Li production. GLAST could detect 
the pionic 7-ray signature from diffuse Galactic emission 
as well as in the EGRB; this would remove the need 
to estimate these components. GLAST should also de- 
tect several Local Group galaxies and thus allow for new 
determinations of the cq smic-ray/star-formation scaling 
llPavlidou fc Fieldsll20r)l . Also,' our calculation is ham- 
pered by lack of evidence of the "pion bump" in the 
Milky Way 7-ray spectrum. Fortunately, it was recently 
demonstrated that future GeV-TeV-PeV gamma-ray ob- 
servations of the diffuse emission from the Galactic Plane 
can determi ne the level of pionic 7 -ray emission in the 
Milky Way l|Prodanovic et alJl20o4 - 

In closing, we have underscored the increasing crucial 
role that 6 Li plays in particle astrophysics and cosmol- 
ogy; the excess we find in solar 6 Li has implications 
throughout these fields. Fortunately, upcoming measure- 
ments of 6 Li and gamma-rays should be able to address 
the questions posed here. 
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